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In class, we saw the connection between the so-called Hermite matrix and the semi- circular law. There 
is actually a deeper story that connects the classical random matrix ensembles to the classical orthogonal 
polynomials studied in classical texts such as [1] and more recent monographs such as [2]. Wc illuminate 
part of this story here. The website www . mathworld . com is an excellent reference for these polynomials and 
will prove handy when completing the exercises. 

In any computational explorations, see if you can spot the interesting feature in the eigenvectors (either 
the first or last row/column) of the corresponding tridiagonal matrix. 



1 Roots of orthogonal polynomials 



Any weight function w{x) on an interval [a,b] (possibly infinite) defines a system of orthogonal polynomials 
7r„, n = 0, 1, 2, . . . , satisfying 



T^mT^nWix) dx = Smn- 

The polynomials can be generated easily, because they satisfy a three-term recurrence 
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(Note that 6o is taken to be zero.) 

Perhaps surprisingly, the three-term recurrence can be used to find the roots of 7r„ as well. Simply form 
the symmetric tridiagonal matrix 

ai bi 
bi a2 62 

62 as 63 



The roots of 7r„ are the eigenvalues of T„! 



bn-2 bn-l 
bn-l CLn 



Remark 1. The symmetry ofT follows from two choices which we made. First, in (1), 7r„ was normalized to 
have unit length under the inner product (/, g) = fgw. Second, (3) was arranged in the form XTTn — RHS. 
Not all authors follow these conventions. 



To see that the eigenvalues are the roots, consider vectors of the form 

V = (7ro(x),7ri(a;), . . . , 7r„_i(a;))'^. 

The recurrence (3) imphes that {Tv)i = XTTi-i{x) = xvi for i = 1, 1, regardless of the value of x. 
Now, if X is a root of TTn{x), then (3) specializes to 

XTTn-iix) = 6„_i7r„_2(a;) + a„7r„_ia; 

i.e. 

XVn = 6„_iW„_i + a„u„ = {Tv)n- 
Therefore, every root of 7r„(a;) is an eigenvalue of Tn- 

2 Hermite polynomials 

Let w{x) = on M. According to mathworld.wolfram.com, the system of polynomials _ff„ defined by 

/oo 
H,nHnw{x) dx = (5,„„2"n!V7r 
-oo 

satisfies the recurrence 



Hn+i{x) = 2xH„{x) - 2nH„-i{x). 
Let the orthonormal Hermite polynomials be defined by Hn{x) = 2nf^J~\. i/i -^^i-'^)' that 



(4) 
(5) 



and substitute 2'^/^{nl)^/^TT^/^Hn{x) for H„{x) in (5) to find 

xHn{x) = -^^/nHn^i{x) + -^y/n + lH„+ii 



V2' 



These equations are exactly of the form (l)-(3), so the eigenvalues of 
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arc the roots of iJ„. 

Exericse: Write a MATLAB function that generates this tridiagonal matrix. Use the eig command to 
compute its eigenvalues for any choice of n . Histogram the roots of _ff„ for n = 50, 100, 250, 300, 500 
using the histn function. Compare it to the semi-circular law. Does it make a difference whether the the 
generalized ensemble is real or complex? 



3 Laguerre polynomials 



Let Wa^x) = a;"e ^ on [0,oo). 

Exercise: Verify that the Laguerre polynomials defined by 



LlKw dx 



satisfy the recurrence 

L^{x) = 



v/r(a + i) 



xL^{x) = -V^\/Q' + 7iL"_i(a;) + (a + 1 + 2n)L^{x) - Vn + iVa + n + 1L^+i(.t). 
Therefore, the roots of are given by the eigenvalues of 



a + 1 -\/TVa + 1 
-V^VaT2 a + 5 



1 Q - 1 + 2n 



Exericse: Write a MATLAB function that generates this tridiagonal matrix. Use the eig command to 
compute its eigenvalues for any choice of n and a. Histogram the roots of L„ for different values of n and 
m using the histn function. Compare it to the Marcenko-Pastur distribution for the generalized Laguerre 
Ensemble. Does it make a difference whether the the generalized ensemble is real or complex? Note: 
a = /3(n — m + l)/2— 1 where /3 = 1 when the elements are real and (3 = 2 when the elements are complex. 
The parameters m and n are associated with the rows and columns of the Gaussian matrix used to 
construct the Wishart matrix. Recall that m/n = c. 



4 Jacobi polynomials 

Jacobi polynomials arc orthogonal with respect to the weight w{x) = (1 — + x)^ on the interval 

Exercise: Verify that the tridiagonal matrix describing the recurrence is produced by the following MATLAB 
code. 



function t = trijacobi(n,a,b) 

7oTRIJACDBI Returns the tridiagonal matrix of polynomial recurrence coefficients 
7„ TRIJACOBI (N, A,B) returns a tridiagonal matrix whose entries correspond 
7„ to the coefficients describing the recurrence between 

7„ the Jacobi Polynomials . 

7. 

7o The eigenvalues of this matrix correspond to the roots of the N-th 

7o Jacobi polynomial with parameters A and B. 

7. 

7o N, A and B are integers. 

7. 
7. 
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7o [2] Alan Edelman, Random Matrix Eigenvalues. 

7o [3] Gabor Szego, Orthogonal Polynomials, American Mathematical 

7o Society, Providence, 1975. 4th Edition. 

7o [4] Eric Weisstein, Jacobi Polynomial." From MathWorld — 

7o A Wolfram Web Resource. 

7o http : //mathworld . wolfram . com/ JacobiPolynomial . htm 

7. 

7. Brian Sutton, Sept. 2004. 

7. $Revision: 1.0 $ $Date: 2004/09/23 11:35:18 $ 
i = (l:n-l)'; 

dl = 2*sqrt(i.*(a+i) .*(b+i) .*(a+b+i) ./(-l+a+b+2*i) ./(l+a+b+2*i)) ./(a+ ... 

b+2*i) ; 

dO = -(a-2-b-2) ./(a+b+2*(0:n-l) O ./(2+a+b+2*(0:n-l) O ; 
t = spdiagsC [[dl;0] dO [0;dl]],[-l l],n,n); 

Code 1 



Exericse: Use the eig command to compute its eigenvalues for any choice of n. Histogram the roots of T„ 
for different values of n, a and b using the histn function. Compare it to the MANOVA density for the 
MANOVA matrix. Hint: You might have to shift the spectrum for it to line up exactly. Does it make a 
difference whether the the generalized ensemble is real or complex? Note: a = fi (ni — to + l)/2 — 1 and 
b = (3 (n2 — TO + l)/2 — 1 where (3 = 1 when the elements are real and (3 = 2 when the elements are 
complex. The parameters to, ni and n2 are associated with the rows and columns of the Gaussian matrices 
used to construct the MANOVA matrix. Recall that m/ni ~ c\ and m/n2 ~ C2. 
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